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ABSTRACT 

We have undertaken a deep, wide-field H I imaging survey of M31, reaching a maximum 
resolution of about 50 pc and 2 km s~^ across a 95x48 kpc region. The H I mass and brightness 
sensitivity at 100 pc resolution for a 25 km s^^ wide spectral feature is 1500 Mq and 0.28 K. Our 
study reveals ubiquitous H I self-opacity features, discernible in the first instance as filamentary 
local minima in images of the peak H I brightness temperature. Local minima are organized 
into complexes of more than kpc length and are particularly associated with the leading edge of 
spiral arm features. Just as in the Galaxy, there is only patchy correspondence of self-opaque 
features with CO(l-O) emission. We have produced images of the best-fit physical parameters; 
spin temperature, opacity-corrected column density and non-thermal velocity dispersion, for the 
brightest spectral feature along each line-of-sight in the M31 disk. Spectroscopically opaque 
atomic gas is organized into filamentary complexes and isolated clouds down to 100 pc. Localized 
opacity corrections to the column density exceed an order of magnitude in many cases and add 
globally to a 30% increase in the atomic gas mass over that inferred from the integrated brightness 
under the usual assumption of negligible self-opacity. Opaque atomic gas first increases from 20 
to 60 K in spin temperature with radius to 12 kpc but then declines again to 20 K beyond 25 kpc. 
We have extended the resolved star formation law down to physical scales more than an order 
of magnitude smaller in area and mass than has been possible previously. The relation between 
total-gas-mass- and star-formation-rate-density is significantly tighter than that with molecular- 
mass and is fully consistent in both slope and normalization with the power law index of 1.56 
found in the molecule-dominated disk of M51 at 500 pc resolution. Below a gas- mass-density of 
about 5 M0pc~^, there is a down-turn in star-formation-rate-density which may represent a real 
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local threshold for massive star formation at a cloud mass of about 5T0^ M, 



Subject headings: galaxies: individual (M31) - galaxies: Local Group - galaxies: ISM - radio lines: 
galaxies - stars: formation ^ 



1. Introduction 

As the nearest major external galaxy and pos- 
sibly the dominant member of the Local Group, 
Messier 31 has received a great deal of atten- 
tion since its' first documented te lescopic obser - 
vation by Simon Marius in 1612 (iMariud ITili ) . 
The many stellar and gaseous constituents of M31 
continue to be studied with ever greater sensitiv- 
ity and resolution, permitting insights that are in 
many cases difficult or even impossible to achieve 
for our own Galaxy given our vantage point within 
the disk. One of these areas, where an external 
yet nearby vantage point is crucial to achieving a 
global overview, is the study of the neutral inter- 
stellar medium. Although we can achieve remark- 
able physical sensitivity and resolution of the neu- 
tral medium of the Galaxy, resolving individual 
structures of lOO's of AU in size and as little as 



several Jupiter masses (|Braun fc Kanekaii 12005.) . 
it is very difficult to estimate even the total mass 
of the neutral medium, let alone its overall distri- 
bution and relationship with other components. 
Studies of the neutral medium of M31 began with 
the pioneering work of I Van de Hulst et al.l (|l957l ) 
as one of the first programs undertaken with the 
Dwingeloo 25m telescope. Subsequent studies 
with larger single dishes and then with interfcro- 
metric arrays have provided enhanced resolution 
and brightness sensitivity, although the highest 
resolution studies have often only been directed at 
a portion of t he disk (eg. , Brinks ^ Shane (1984. ), 



BraunI (|l99lf l. ICarignan et all (|2006l )). 



In this paper we describe our Westerbork and 
Green Bank Telescope program to image the en- 
tire atomic disk of M31 with a uniformly high 
sensitivity and resolution. The mo st recen t effort 
on a comparable scale was that of lUnwiiJ ( 1983h 



which achieved arcmin spatial and 16 km ve- 
locity resolution with an RMS noise of 3.6 K in 
two elliptical fields covering about 12 square de- 
grees. By contrast we will present imaging down 
to 15" spatial and 2.3 km s^^ velocity resolu- 
tion (sampled at 2 km s^^) with a uniform 2.7 K 
RMS over a region of about 24 (7x3.5) square 
degrees. This corresponds to a maximum linear 
resolution of about 50 pc over a total extent of 
more than 90 kpc. With even modest spatial 
smoothing we reach sub-Kelvin brightness sensi- 
tivity (for example 0.37 K at 1 arcmin), making 



our database the most detailed documentation yet 
of the neutral medium of any complete galaxy 
disk, including our own. For comparison, the 
Leiden/ Argentine/Bonn Survey of Galactic H I 
has achieved 36 arcmin and 1.3 km s"^ resolu- 
tion over 360 degrees and reaches 0.08 K RMS 



Kalberla et al. 



( 20051) and the La rge Magellanic 
Kim et al. ( 20031 ) achieves 1 ar- 



Cloud survey of 
cmin and 1.7 km s~^ resolution over 12x10 de- 
grees and reaches 2.5 K RMS. 

We begin with a description of the observations 
and their reduction in ^ and continue with an 
overview of the results in ^ In fJHwe present some 
initial analysis of the data, beginning with a dis- 
cussion of H I self-opaque features and followed by 
a detailed fitting of all high signal-to-noise spectra 
to determine the best fitting combination of the 
physical parameters along each line-of-sight: cool 
component temperature, total column density and 
non-thermal velocity dispersion. We continue with 
the analysis of radial dependence for the relevant 
physical parameters and finally we focus on the 
relationship between surface density of gas mass 
and star formation rate. We close with some brief 
conclusions in ^ We assume a dist ance to M31 
of 785 kpc (jMcConnachie et al.|[2005h throughout. 

Further papers in this series will include: (1) a 
detailed kinematic- and mass-model of M31, (2) 
comparative studies of the morphology of the in- 
ner and outer disks and (3) a study of small-scale 
kinematic features such as outflows and other pe- 
culiar velocity components. 

2. Observations and Reduction 

Interferometric observations of the extended 
M31 H I disk with the Westerbork Synthesis Ra- 
dio Telescope (WSRT) array were obtained on a 
Nyquist sampled (15 arcmin spacing) pointing grid 
defined in the M31 (major, minor) axis coordinate 
system. The disk extent was first determined in a 
precursor total power survey which fully sampled 
a 6x6 degree field at low resolution (36' FWHM) 
using the 14 individual telescopes of the WSRT 
array. A total of 163 synthesis pointing positions 
were chosen from the grid to provide good cover- 
age of the region of anticipated H I emission. Ob- 
servations were then obtained for 27 tracks of 12 
hour duration with the WSRT array, each sam- 
pling 6 (or in one case 7) different pointing po- 
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sitions. These data were observed between Au- 
gust 2001 and January 2002. Each 12 hour track 
involved cycUng through the 6 relevant pointing 
positions with a series of 10 minute snap-shots. 
The tracks were preceeded by half hour observa- 
tions of the calibration sources 3C286 and CTD93 
and followed by half hour observations of 3C147 
and CTA21. A total band- width of 5 MHz (cov- 
ering heliocentric radial velocities between —828 
and -1-227 km s~^) was observed with 512 spec- 
tral channels in two linear polarization products 
yielding a channel width of 2.06 km s~^ and res- 
olution of 2.27 km s~^. Calibration and flagging 
of the data were done in Classic AIPS. After ex- 
ternal band-pass and gain calibration, each point- 
ing was self-calibrated using the continuum emis- 
sion which happened to be present. The contin- 
uum model was subtracted from the visibilities 
and residual continuum was removed with UVLIN. 

Sensitive total power data covering a 7x7 de- 
gree field was acquired with the Green Bank 
Telescope (GBT) telescope over six nights in 
September 2002. These data, with 9.'1 FWHM 
X 1.03 km s~^ resolution and an RMS sensitivity 
in H I column density over relevant linewidths of 
about 5x10^^ cm~^, have already been presented 



Thilker et all (|2004[ ) 



A joint deconvolution of all 163 pointing s was 



done with the Miriad ( Sault et al. I Il995f ) task 



MOSSDI, after forming a weighted sum of the lin- 
early combined synthesis image (employing uni- 
form weighting followed by a Gaussian taper) with 
the relevant GBT total power image. Each of 
the 163 synthesized beams was also modified by 
addition of a Gaussian approximation of the to- 
tal power beam using the same relative weight- 
ings as used when summing the synthesis and to- 
tal power images. The relative weightings were 
determined empirically by minimizing the "short- 
spacing bowl" effect in the synthesis beams, but 
they correspond roughly to the inverse of the rela- 
tive beam areas. The exact choice of weighting is 
not critical, in theory, since the process of decon- 
volution should yield an "ideal" final PSF if all 
emission has been successfully incorporated into 
the component model. In practise it is never pos- 
sible to incorporate all emission into the compo- 
nent model, so it is important to make the "dirty" 
PSF as clean as possible (particularly with re- 
spect to diffuse emission) by an optimized choice 



of weights. 

The linearly combined "dirty" mosaic images 
were multiplied with the combined WSRT sur- 
vey sensitivity pattern in order to weight down 
the noise around the periphery of the sampled re- 
gion before deconvolution. If this were not done, 
then the noisy periphery would be severely over- 
CLEANed relative to the central regions, yielding 
greatly reduced image fidelity. The assumption 
underlying this spatial tapering is that the region 
of significant emission is contained within the sam- 
pling pattern. We will comment below on the ex- 
tent to which this assumption was met in practise. 
The 300 semi-independent spectral channels (each 
of 2 km s^"'^ width) containing M31's H I emis- 
sion (heliocentric velocities of about —621 to —20 
km s~^) were processed individually. Foreground 
H I emission from the Galaxy is present in the field 
of M31 at velocities between about —130 and -1-45 
km s-^ Reduction of the WSRT data over the full 
velocity range of Galactic emission was not under- 
taken at this time. The GBT total power images 
at velocities between —130 and —20 km s~^ were 
modified by subtraction of the median value deter- 
mined outside of the region occupied by the M31 
disk before combination with the WSRT mosaic 
to simplify the deconvolution problem. A positive 
pedestal of emission underlying the entire galaxy 
is a needless complication that undermines image 
fidelity in the deconvolved result. 

Given the non-linear nature of the deconvo- 
lution process and the desire to probe a vari- 
ety of different output spatial ("fuir'=18'.'4xl4'.'8 
(NSxEW), 30", 60"and 120") and velocity (2.3, 
6.2, 10.3, and 18.5 km s^^) resolutions, we carried 
out the joint deconvolution independently for each 
different combination of these parameters. De- 
convolution was in all cases done to a depth of 
twice the theoretical RMS noise level at each res- 
olution. A bug in the MOSSDI source code was 
remedied (affording it protection against a neg- 
ative loop gain). A small adjustment was also 
made to the primary beam model employed within 
Miriad for the WSRT telescopes. It was found em- 
pirically, that a small increase (by 5 %) in the 
assumed FWHM of the analytic approximation 
{cos^{C' ■ V ■ r) for a constant C at frequency v 
in terms of the radius r) for the beam pattern 
yielded the highest dynamic range in the mosaic 
result. The slightly larger than expected beam size 
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(2110 arcsec FWHM rather than 2013 arcsec at 
1420 MHz ) has since been verified by direct mea- 
surement ( Popping &: BraunI 2008h . The resuhing 
RMS noise levels at various spatial and velocity 
resolutions are listed in Tabled] Note that unlike 
the case of total power observations of Galactic 
H I, the detected M31 H I emission in our synthe- 
sis images does not contribute significantly to the 
measured antenna temperature of each dish, and 
therefore does not give rise to an increase of the 
system temperature. 

A locally modified version of Miriad was em- 
ployed that contains many enhancements for high 
performance computing. It extensively exploits 
parallelism at several levels: multi-threading (on 
machines that have multiple CPUs) and vector- 
ization (using SSE (Streaming SIMD Extensions) 
instructions which are available on most contem- 
porary processors). Furthermore, we used task- 
level parallelism by dividing different reduction 
pipelines (specific resolutions and spectral chan- 
nels) over multiple machines. The "Netherlands 
Grid" infrastructure provided us with the re- 
sources to do the computational work. Using these 
techniques, we were able to perform all reductions 
(16 different resolutions, 300 semi-independent 
spectral channels) in 36 hours of elapsed time. A 
previous attempt to carry out the same process- 
ing, before code optimization and using only a 
single high-end server, required several months of 
elapsed time. 

Although data products at all 16 different res- 
olutions were generated directly from the visibil- 
ity data, we also produced velocity-averaged data 
products at a fixed spatial resolution by averag- 
ing the individually processed channels after de- 
convolution. Superior dynamic range was found 
for these post-deconvolution averaged products at 
spatial resolution of both 60 and 120 arcseconds, 
while the directly imaged products were supe- 
rior at full and 30 arcsecond resolutions. Ap- 
parently, different deconvolution limitations dom- 
inate in these two regimes. At high spatial resolu- 
tions, the increased signal-to- noise ratio afforded 
by visibility-based velocity smoothing improves 
the recovery of extended emission structures, while 
at low spatial resolutions we already encounter the 
finite image fidelity of the deconvolution in an in- 
dividual velocity channel. In this case, averaging 
spectral images after deconvolution can improve 



the dynamic range of the result. The nominal 
RMS noise levels after spectral smoothing scale 
simply with the square root of the eventual veloc- 
ity resolution. H I column density sensitivity at an 
appropriate spectral smoothing for faint emission 
features (of 18.5 km s~^) is tabulated in Table [1] 
at a series of spatial resolutions. 

Images of the peak observed line brightness and 
the integrated brightness and associated mean ve- 
locity (the 0*'* and 1^* moments) were produced at 
each resolution after first masking the data-cube. 
The mask was constructed from the cube of 120 
arcsec and 18 km s^^ resolution, by first blank- 
ing brightnesses below 1.5cr and then interactively 
masking with an arbitrary polygon the "island" 
of most significant emission at each velocity. In- 
teractive isolation of M31 H I emission was most 
critical in the velocity range that overlaps that of 
the Galaxy (—130 to —20 km s^^). Subtraction 
of the off-M31 median total power before decon- 
volution already served to place the M31 emission 
on a plausible base level. Only a single contiguous 
(both spatially and in velocity) "island" of emis- 
sion at each velocity was included in the blanking 
mask. Although isolated features in the field may 
be associated with M31, we are unaware of how 
to distinguish them from features in the Galactic 
foreground. All features in the field which were not 
contiguously connected with the M31 disk emis- 
sion were presumed to be part of the Galactic fore- 
ground and were excluded. 

Images of integrated brightness were rescaled 
to represent column density assuming a negligible 
H I line opacity. Although this is the usual as- 
sumption, it is clear that it is incorrect in practise. 
Braun fc Walterbod (|l992l) have applied a simple, 
single temperature opacity correction to the ob- 
served H I brightnesses from the North-East half 
of M31 and find local corrections of a factor of 
2 and an integrated correction of some 20%. We 
address the issue of H I opacity in greater detail 
below by directly fitting physical models to indi- 
vidual spectra for all lines-of-sight. 

Another useful type of image for diagnostic pur- 
poses is that o f "velocity coh erence" (VC) (as first 
introduced by iBrauii ( 1997 )). defined as the ratio 
of peak brightness at two different spectral reso- 
lutions. For example, YC^Q/^lx^y) = P3o{x,y) 
j P[^x,y\ where Pii{x,y) is the peak brightness 
observed at a spectral resolution of R km s^^. In 
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Table 1 

RMS Sensitivity at Various Resolutions 



Resolution RMS Sensitivity 



Beam 


Ay 


AS 


ATb 


ANhi 


AMhi 


[arcsec] 


[km s-i] 


[mJy/Bcam] 


[K] 


[10i8cm-2] 


[Mq] 


18x15 


2.3 


1.2 


2.7 


11.3 


390 


30 


2.3 


1.5 


1.0 


4.2 


490 


60 


2.3 


2.2 


0.37 


1.6 


720 


120 


2.3 


4.2 


0.15 


0.63 


1400 


18x15 


18.5 


0.42 


0.95 


32 


1100 


30 


18.5 


0.53 


0.35 


11.8 


1400 


60 


18.5 


0.78 


0.13 


4.4 


2100 


120 


18.5 


1.5 


0.053 


1.8 


4000 



the absence of noise, < VC < 1, when compar- 
ing coarser with finer spectral resolutions, since 
an infinitely broad spectral feature will always re- 
tain the same peak brightness (VC = 1), while an 
extremely narrow feature will be strongly diluted 
(VC 0). Images of velocity coherence are a ro- 
bust measure of linewidth, since no assumption 
regarding the intrinsic line-shape, nor reliance on 
a fitting algorithm are required for it's calcula- 
tion. Moreover, even when spectra become multi- 
valued, as they do in highly inclined (like M31) or 
kinematically disturbed systems, the VC can con- 
tinue to probe the attributes of the dominant spec- 
tral feature along each line-of-sight rather than a 
smeared version of all spectral content. In the par- 
ticular case of an intrinsically Gaussian line pro- 
file with dispersion, ct, that is fully resolved at the 
higher resolution and then convolved with a Heav- 
iside (top-hat) function of width, AV, one gets the 
relation, a = AV ■ VC/[{1 - VC)V2^]. 



3. Overview of Results 



An overview of the joint dcconvolution result is 
shown in Figure [1] where the WSRT mosaic sensi- 
tivity pattern is overlaid on a low resolution (120 
arcsec, 10 km s~^) i mage of peak H I brightness. 
The Karma package ( Gooch 19961 ) has been em- 
ployed for many of the illustrations. Sensitivity 
contours are shown at 10 and 50% of the peak 
sensitivity which results from the sampling pat- 



tern. The combined sensitivity pattern has what 
turn out to be unfortunate local minima on cither 
side of the minor axis. Although little evidence for 
emission was discerned at these locations in our 
pilot WSRT total power observations, it is clearly 
present in the more sensi tive total pow e r data that 
we have acquired since (jThilker et al.l (|2004l ) and 
Braun fc Thilked (j2004l )). Diffuse patches of emis- 



sion are apparent on both sides of the minor axis. 
The apparent termination of these minor axis ex- 
tensions is likely influenced by the multiplication 
of the brightness distribution with the sensitivity 
pattern (to control noise amplification during de- 
convolution, as outlined in fQ] above). More ex- 
tensive sampling of the minor axis with additional 
mosaic pointings will be required for high fidelity 
imaging at arcmin resolution. Our total power 
data (with up to 9 arcmin resolution) do not suf- 
fer from this limita tion. Similar a ttenuation af 
fects Davies' Cloud (|Daviesl (|l975l ). IPe Heii et al 
(|2002l )) the diffuse structure at 5=42? 5 to the 
northwest of the M31 disk. The southernmost of 
the H I clump s at (5=38?8 and due South of M31 



([Thilker et al.ll2004l ) is also somewhat attenuated. 
Sensitive poi nted observations of the se features are 
presented in lWestmeier et al.l (|2005l ). 

Greater detail can be discerned in Figure [Din 
which peak brightness is shown as determined at 
60 arcsec (220 pc) and 6 km s~^ resolution, al- 
though even at this modest resolution, the PSF 
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Right Ascension (J2000) 

Fig. 1. — Peak brightness of H I emission determined at 120 arcsec and 10 km s~^ resolution is overlaid with 
the combined sensitivity pattern of the 163 pointings in the mosaic observation. Contours are drawn at 10 
and 50% of the peak sensitivity. Peak brightness is shown on a logarithmic scale which saturates at 60 K. 
The beam FWHM is indicated in the lower left corner. 
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(plotted in the lower left corner of the image) is 
barely resolved on a printed page, given that the 
image spans 111 kpc across the diagonal. Zoom- 
ing in at this resolution on the central disk as in 
Figure [3] finally permits appreciation of which fea- 
tures are resolved. Jumping now to 30" (110 pc) 
resolution on the same field in Figure |4] begins 
to reveal the intricate filamentary structure of 
the neutral ISM. Point-like local minima are seen 
toward the most prominent background sources 
(J004218-f412926 and J004648-H420855) against 
which H I absorption is detected. Increasing the 
angular resolution to the maximum available of 
about 55 pc (18.4x14.8 arcsec) as in Figure [5] 
does not appear to appreciably add detail to the 
filamentary structure, while increasing the noise 
level significa ntly. This is consistent with what 
was found by iBraun ( 1995 ) when analyzing H I 
emission data from the northeast half of the in- 
ner M31 disk; beginning with data having 10" 
and 5 km s~^ resolution, spatial smoothing to 40" 
resolution was required before the mean peak de- 
tected brightness decreased by 8%, while velocity 
smoothing to even 10 km resulted in a mean 
decrease of peak brightness of 15%. The major H I 
cloud complexes which are discernable in external 
galaxies with current instrumentation are essen- 
tially resolved with about 100 pc and 5 km s~^ 
resolution. However, it must also be stressed that 
significant structure is present at even the highest 
angular and velocity resolution that is available. 
We will briefly discuss some examples in 21 below. 

Continuing on from the best signal-to-noise in- 
ner disk overview of Figure IH further magnifica- 
tion is required to discern individual features on a 
printed page. The three images in Figs. [MSI give 
some appreciation of the wealth of detail which 
can be discerned in the inner disk of M31. Peak 
H I brightness is shown in the left hand panels 
of each figure, while velocity coherence VCig. 5/2.1 
(defined in Jj2]) for the same region is shown on 
the right. A striking aspect of these images is the 
appearance of filament-like minima as well as posi- 
tive morphological features, as seen in peak bright- 
ness. While positive emission features are not sur- 
prising when imaging the H I distribution of an 
external galaxy, the filamentary minima are quite 
novel. A further striking feature is the correspon- 
dence of high velocity coherence with the majority 
of these local minima. Apparently the filamentary 



features having depressed peak brightness are of- 
ten accompanied by a substantial broadening of 
the underlying spectral profile. This trend is il- 
lustrated in Figs. HI and [TO] (at 30"and 2.3km s^^ 
resolution), where sequences of spectra are shown 
which contrast the appearance of these local min- 
ima (central panels) with flanking spectra which 
are offset perpendicular to each "dark" filament by 
only about 30 arcsec (left and right panels). The 
locations of these cross-cuts are shown in Figs. [Bl 
[7] and [8] From these spectra it is apparent that 
the line profiles toward such local minima in peak 
brightness are not only broad, but have a distinc- 
tive fiat-topped appearance that is strongly sug- 
gestive of self-opacity in the emission line. The 
smooth curves overlaid on each spectrum are the 
parameterized model-fits to each spectrum that we 
describe in detail below in §4.21 Note that the 
RMS sensitivity of these spectra, is that given in 
Table [T] of about 1 K, implying that some of the 
fluctuations seen toward the flat-topped peaks are 
quite significant. The parameterized fits are in- 
deed consistent with a significant self-opacity of 
the H I emission within this entire system of fila- 
mentary local minima. 

4. Discussion 

4.1. Discrete Opaque Features 

The H I self-absorption (or HISA) phenomenon 
has been well-doc umented in the Galaxy since the 
pioneering work of RadhakrshnanI (119601 ) and more 
recently by iGibson et all (|2005l ) in a large swath 
of the Galactic plane (145° > I > 75°). The 
physical conditions required to witness HISA are a 
"sandwich" geometry in which a cooler spin tem- 
perature foreground gas is seen toward a higher 
brightness temperature background gas at the 
same radial velocity. The brightness decrement of 
a HISA feature is proportional to this temperature 
difference in the absence of additional complicat- 
ing factors. The simplest circumstance to imagine 
is a single warm, semi-opaque background feature 
which has a brightness temperature comparable 
to it's own spin temperature together with a sin- 
gle cooler, opaque feature in the foreground. In 
this simple case the HISA temperature decrement 
would give some information about the spin tem- 
perature of the foreground component. More com- 
plicated geometries are easily conceivable. Adding 
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Right Ascension (J2000) 



Fig. 2. — Peak brightness of H I emission determined at 60 arcsec and 6 km s^^ resolution. Peak brightness 
is shown on a logarithmic scale which saturates at 90 K. The beam FWHM is indicated in the lower left 
corner. 
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Right Ascension (J2000) 

Fig. 3. — Peak brightness of H I emission determined at 60 arcsec and 6 km s^^ resolution of the central 
50% of the survey region. Peak brightness is shown on a square-root scale which saturates at 90 K. The 
beam FWHM is indicated in the lower left corner. 
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Right Ascension (J2000) 

Fig. 4. — Peak brightness of H I emission determined at 30 arcsec and 6 km s^^ resolution of the central 
50% of the survey region. Peak brightness is shown on a square-root scale which saturates at 120 K. The 
beam FWHM is indicated in the lower left corner. 
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Fig. 5. — Peak brightness of H I emission determined at about 15 arcsec and 6 km s~^ resolution of the 
central 50% of the survey region. Peak brightness is shown on a square-root scale which saturates at 120 K. 
The beam FWHM is indicated in the lower left corner. 
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Right Ascension (JaOOO) Right Ascension (jaOOO) 



Fig. 6. — Peak brightness (left) and velocity coherence (right) of H I emission determined at 30 arcsec 
resolution of the North-East disk region. Peak brightness is shown on a square-root scale which saturates at 
120 K. Velocity coherence (defined in the text) is measured between 18.5 and 2.1 km s~^ spectral smoothing. 
The black bars indicate the location of the spectra shown in Fig. [9l The beam FWHM is indicated in the 
lower left corner. 
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Right Ascension (J^OOO) Right Ascension (jaOOO) 

Fig. 7. — Peak brightness (left) and velocity coherence (right) of H I emission determined at 30 arcsec 
resolution of the central survey region. Peak brightness is shown on a square-root scale which saturates at 
120 K. Velocity coherence (defined in the text) is measured between 18.5 and 2.1 km s~^ spectral smoothing. 
The black bar indicates the location of the spectra shown in Fig. [9l The beam FWHM is indicated in the 
lower left corner. 
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Right Ascension (JaOOO) Right Ascension (jaOOO) 



Fig. 8. — Peak brightness (left) and velocity coherence (right) of H I emission determined at 30 arcsec 
resolution of the South- West disk region. Peak brightness is shown on a square-root scale which saturates at 
120 K. Velocity coherence (defined in the text) is measured between 18.5 and 2.1 km s~^ spectral smoothing. 
The black bars indicate the location of the spectra shown in Fig. [TOl The beam FWHM is indicated in the 
lower left corner. 
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Fig. 9. — H I spectra of dark filaments. The three sequences of spectra (top, middle and bottom) each repre- 
sent a cross-cut perpendicular to a filamentary local minimum in the image of peak brightness temperature 
shown in Figs. [6]and[7l The spectrum in each central panel is in the direction of the dark filament, while 
the left and right spectra are offset by about 30 arcsec (as indicated by the coordinates). The resolution is 
30"and 2.3 km s^^. Note the depressed, fiat-topped profiles toward the filament, flanked by brighter, more 
smoothly peaked profiles. The smooth curve overlaid on each spectrum is the parameterized fit. 
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Fig. 10. — H I spectra of dark filaments. The three sequences of spectra (top, middle and bottom) each 
represent a cross-cut perpendicular to a filamentary local minimum in the image of peak brightness temper- 
ature shown in Fig. [H The spectrum in each central panel is in the direction of the dark filament, while 
the left and right spectra are offset by about 30 arcsec (as indicated by the coordinates). The resolution is 
30"and 2.3 km s^^. Note the depressed, fiat-topped profiles toward the filament, fianked by brighter, more 
smoothly peaked profiles. The smooth curve overlaid on each spectrum is the parameterized fit. 
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Fig. 11. — Peak brightness of H I emission determined at 30 arcsec and 6 km s~^ resolution of the North- 
East disk region with contours of integrated CO(l-O) (from (i Nieten et al. 2006h ) superposed at 2.34 x 1.2" 
K-km/s, for n=l,2,3. . .. 
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Fig. 12. — Peak brightness of H I emission determined at 30 arcsec and 6 km s ~^ resolution of the South- 
West disk region with contours of integrated CO(l-O) (from (i Nieten et al. 20061 )) superposed at 2.34 x 1.2" 
K-km/s, for n=l,2,3. . .. 



18 



a third warm, semi-opaque foreground component 
to make a warm-cool-warm "sandwich" would sig- 
nificantly fill-in the detected temperature decre- 
ment. Clearly, the interpretation of HISA tem- 
perature decrements in physical terms is poorly 
constrained. 

The large st individual features seen in the 



Gibson et al.l (|2005l ) Galactic HISA sample ex- 
tend over about 5° at 2 kpc distance, and so 
have lengths as large as 175 pc. A subset of the 
Galactic HISA features, such as those kinemati- 
cally associated with the Perseus spiral arm, are 
organized over tens of degrees making the com- 
plexes at least 1 kpc long. The departure from 
the Galactic ed ge-on geomet ry to the 78° incli- 
nation of M31 (|Braunlll991l) means that the nec- 
essary alignment conditions for witnessing HISA 
are essentially eliminated on large scales in M31. 
Instead, the "sandwich" geometry that can yield 
large-scale HISA features in the Galaxy would be 
projected into spatially resolved, parallel "slices" 
of semi-opaque gas of different spin temperature. 
We identify the filamentary local minima in peak 
brightness seen in M31 with the colder opaque 
features that are responsible for large-scale HISA 
in the Galaxy. We use the term "self-opaque" to 
emphasis the importance of internal optical depth 
effects in determining the profile shapes of these 
features, as distinct from "self-absorption" which 
implies a substantial temperature contrast of fea- 
tures that overlap both along the line-of-sight and 
in radial velocity. 

By comparison to the kpc extent of Galactic 
HISA complexes, the self-opaque filamentary min- 
ima in M31 (for example the one running from 
(a, (5) - (00:45:25, +41:36) to (00:45:55,+41:50) 
in Fig [6]) are often in excess of 10 arcmin in 
length corresponding to more than 2 kpc. One 
such linear feature is seen to cross very near the 
line-of-sight to the background continuum source 
J004218-f 412926 (=B0039+412), as seen in FigEl 
H I absorption measurements for this source and 
several others have been published previously in 
Braun fc WalterbosI ( 1992 ) . We will comment fur- 
ther on these sources below. More complex lo- 
cal minima of all sizes can be seen in the M31 
data wherever sufficient local intensity is present 
in both position and velocity. In contrast, the fila- 
mentary network of systematically broadened pro- 
files (traced by the velocity coherence) are even 



more ubiquitous, since profile broadening can be 
detected independent of the emission characteris- 
tics of the surroundings. 

A subset of Galactic HISA features is closely 
associated with molecul ar clouds as trace d by O H 
and CO emission (eg. IColdsmith fc Lil (|2005l )). 
while th e majority have no apparent CO coun- 
terparts ( Gibson et al.l 2000 ). The co njecture has 
been made (eg. Gibson et al. ( 20051 )') that those 
HISA features without associated CO detections 
may represent an earlier evolutionary phase in 
molecular cloud formation and that such features 
may be organized over kpc-scale regions by the 
passage of spiral density wave (or other forms 
of) shocks. Only after subsequent cooling and 
compression would complex molecule- and star- 
formation proceed. Some of the most dramatic 
self-opaque features discernible in the M31 peak 
brightness images appear to lie on the leading 
edges of spiral arm structures, for example the fea- 
ture extending from {a,S) -(00:45:25,-1-41:36) to 
(00:45:55,-1-41:50). (We can deduce the upstream 
edge of bright H I features from the counter- 
clockwise sense of rotation that follows from the 
approaching major axis PA ~ 232° East of North 
and the near-side minor axis PA = 322°). This 
appears to be quite a general phenomenon; dis- 
cernible as a sharp local minimum in peak H I 
along the leading edge of essentially every arm seg- 
ment with a suitably bright background. Compar- 
ison of the peak H I brightness with the integrated 
CO(l-O) obser ved with the IRAM 30m at 23 arc- 
sec resolution (jNieten et al.ll2006l) in Figs. [TTHT^ 
suggests a rather good correlation of the peak 
brightness of H I with integrated CO at this res- 
olution of about 100 pc at radii less than about 
12 kpc; with a distinct decline of associated CO 
to larger radii. In contrast, the filamentary self- 
opaque minima, particularly those associated with 
the leading edges of spiral features, have only 
very patchy correspondence with CO(l-O) fea- 
tures, even near 12 kpc. 

4.2. Fitting for Physical Parameters 

While some indications for self-opacity in the 
H I line are indicated by qualitative features, 
like flat-topped spectral profiles, obtaining quan- 
titative estimates of the "hidden" mass in neu- 
tral hydrogen is considerably more challenging. 
The "classical" method of H I opacity correc- 
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tion (ISchmidd (Il957l). iHenderson et alj (|l982f ). 
iBraun fc WalterbosI tmH ) has consisted of cal- 
culating the column density from: 



Such a n approach h a s been considered previ- 



ously by iRohlfs et al 



HI 



-1.823 X 10^^TsAVln{ 1 - Icm"^ 



(1) 

where a single cool component spin temperature, 
Ts, is assumed to apply for all spectra. The rele- 
vant spin temperature has been chosen to be sim- 
ilar to the highest brightness temperature actu- 
ally obs erved, for example = 125 K for the 
Galaxy (jHenderson et al.lll982n . or has been esti- 
mated from a fit to the relation between observed 
absorption opacity an d the emission brightness , 
Ts = 175 K for M31 (|Braun fc WalterbosI 1 199^) . 
In both cases, a single spin temperature is assumed 
to apply throughout the galaxy in question, which 
is clearly a poor assumption, and the highest ob- 
served brightnesses (those with Tg > Ts) have an 
undefined column density. 

Rather than assuming that a single spin tem- 
perature might be appropriate for an entire galaxy, 
let us consider a more realistic scenario in which 
the major atomic clouds in a galactic disk are ap- 
proximately isothermal on scales of 100 pc. From 



BraunI (|1997f ) we have an estimate of the face-on 



surface covering factor of such features within the 
star forming disk in a sample of nearby galaxies of 
about 15%. This implies that even at the moder- 
ately high inclination of M31, of 78°, we are un- 
likely to have more than a single major cloud along 
any line-of-sight in the absence of strong warping. 
(Even in the strongly warped case, there need not 
be widespread velocity overlap of these compo- 
nents.) In addition, the good correlation between 
absorption opacity and the emission brightness 
iBraun fc Walterbos 1992j in both the Galaxy and 
M31 suggest that the contribution of warm (5000 - 
10"*^ K), optically thin gas in the extended environ- 
ment of major clouds to the emission brightness is 
only a few Kelvin, although distributed over the 
relevant line-width of some lO's of km s~^. These 
considerations suggest that detailed modeling of 
the peak (and not the wings) of an emission profile 
as an isothermal cloud may provide a useful esti- 
mate of the underlying physical parameters, pro- 
vided the physical and velocity resolution is suf- 
ficiently high (about 100 pc and 2 km s^"'^) and 
the galaxy is viewed at an inclination of less than 
about 81°. 



(|1972I ). who have docu- 
mented the precision and expected biases that 
might pertain to parameter estimation as a func- 
tion of the peak signal-to-noise for an unblended, 
isothermal cloud. We will return to the question 
of the expected uncertainties in such an approach 
below. 

The line profile due to a spatially resolved 
isothermal H I feature in the presence of turbulent 
broadening can be written as 



Tb{V) = Ts{l - exp[-T{V)]} 



(2) 



with 



5.49 X / 

where the velocity dispersion, cr, has units of 
km s^^ and is the quadratic sum of an assumed 
thermal and nonthermal contribution = (4 + 
'^nt) with the thermal contribution given by ctt = 
0.093-\/T'fe, for a kinetic temperature, Tfe. While 
such profiles have a Gaussian shape for low ratios 
of the H I column density (in units of cm^'^) to spin 
temperature, Nhj/Ts they become increasingly 
flat-topped when this ratio becomes high. For sim- 
plicity we assume that the H I spin and kinetic 
temperatures are equal, Tk = Ts, which should 
apply to regions of mod erate volume density al- 
though not in general, cf. Lisztl ( 2001 ). We do not 
consider the possibility of self-absorption (as dis- 
cussed above in S I4.1|) due to multiple temperature 
components along individual lines-of-sight, despite 
seeing occasional evidence for this phenomenon on 
small scales. Simple simulations to explore such 
multi-temperature models demonstrate the great 
challenges of meaningfully constraining the ab- 
sorber properties. Even negligible masses of cold 
H I can result in dramatic modulation of the com- 
posite spectrum. The additional free parameters 
associated with such a composite spectrum make 
for an intractable fitting problem. 

We show examples of model isothermal (half- 
)spectra in Fig. [T31 Shown in the figure are a se- 
quence of increasing aNT = 2 to 12 by 1 km s~^ at 
a fixed log(NH/) = 22 and log(Ts) = 1.7 (on the 
left) as well as a sequence of increasing log(Ni//) 
= 21.5 to 23.5 by 0.2 at a fixed log(Ts) = 1.7 and 
o'NT = 7 km s^^ (on the right). We found the 
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Fig. 13. — Illustration of model spectra. Shown is a sequence of "half-" spectra with increasing unt = 2 to 
12 by 1 km s~^ at a fixed log(N///) = 22 and log(Ts) — 1.7 (on the left) as well as a sequence of increasing 
log(Nff/) = 21.5 to 23.5 by 0.2 at a fixed log(T5) = 1.7 and a^r = 7 km s'^ (on the right). 



best-fitting model spectra for the brightest spec- 
tral feature along each line-of-sight using the 30" 
data with full velocity resolution when the ob- 
served peak brightness temperature exceeded 10 K 
(10 times the RMS noise) and the truncated peak 
(where brightnesses greater than 50% of the peak 
were encountered) spanned at least five indepen- 
dent velocity channels. This truncation was done 
to both isolate single spectral components from 
possibly blended features as well as eliminating 
potential broad wings from the fit. The data were 
compared to a pre-calculated set of model spectra 
spanning the range log(Njj/) = 20.0 to 23.5 by 
0.01, log(Tc) = 1.2 to 3.2 by 0.01 and \og{aNT) 
= 0.3 to 1.5 by 0.04. A search in velocity offset was 
done over displacements of —4 to -1-4 by 1 km s~^ 
with respect to the line centroid as estimated by 
the first moment of each truncated peak. We pur- 
posely limited the range of log(Tc) to that ex- 
pected for the Cool Neutral Medium (CNM) where 
opacity effects might be expected to occur, but 
included sufficiently high temperatures (at least 
ten times the observed peak brightness tempera- 
ture) that negligible opacity solutions are always 
available for the fit. Lines-of-sight which might be 



dominated by the Warm Neutral Medium (WNM) 
with spin temperature of perhaps 8000 K, will be 
fit by a moderately high T5 (since such profiles 
would not display opacity effects), combined with 
an apparent "non-thermal" dispersion of about 



10 km s 



There is no method to discriminate 



between such intrinsically warm gas and cooler gas 
with a comparable non-thermal broadening of the 
profile. All that can be discriminated from the 
profile shape in the first instance is the presence 
or absence of opacity effects. Only if significant 
opacity effects are detected (implying r > 1) can 
the physical parameters be constrained. However, 
since our primary goal with this procedure is the 
correction of apparent column density, this is pre- 
cisely the regime in which we can expect some util- 
ity of this fitting method. This conclusion is con- 



sistent with the experience of lRohlfs et aTl (|l972h 



who indicate that useful parameter fits are only 
achieved with a peak signal-to-noise greater than 
about 10 (such as we require in our fitting) paired 
with peak opacities in excess of unity. Much higher 
signal-to-noise is required for a meaningful fit at 
lower peak opacities. Systematic biases of the fit 
parameters are expected to be below about 10% 
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Fig. 14. — Illustration of spectral fits. Actual spectra at the positions indicated above each box with 
best- fitting models overlaid (smooth curves). The goodness-of-fit is indicated for each spectrum with the 
normalized . Goodness of fit decreases from left-to-right and top-to-bottom, primarily due to the effects of 
velocity blending of physically unrelated features. A cut-off ^ 25 was found to best distinguish blended 
from unblended spectra. 
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in this regime. 

Plausible spectral fits were possible for most 
lines-of-sight. However, there are locations in 
the M31 disk where emission features which are 
likely at significantly different distances along the 
line-of-sight become strongly blended in veloc- 
ity. Such circumstances are most apparent in 
position-velocity images made parallel to the ma- 
jor axis, in which continuous spiral arm segments 
can be tracked as linear features with approxi- 
mately const ant orientation (as documented by 
BraunI (|l99l[ )). The gradient of velocity with ma- 
jor axis distance. 



dV _ Vc(R)sin{i{R)) 
dX ^ R 



(4) 



is proportional to the local projected rotation ve- 
locity divided by the radial distance of the feature 
within the galaxy. The ocurrence of apparently in- 
tersecting linear features in position- velocity plots 
is a strong indication for velocity blending of phys- 
ically unrelated features which reside at different 
radii in the galaxy. The most prominent instances 
of such velocity overlap occur in four locations 
within M31 at major axis distances of both ±30 
arcmin and minor axis distances of both ±10 ar- 
cmin. 

The challenges of velocity blending are illus- 
trated by the sequence of individual spectra with 
overlaid fits shown in Fig. [TH The sequence of 
spectra from left to right and top to bottom are 
ordered by decreasing quality of fit, as measured 
by the normalized = T,i{Di-M,)'^/{Na'^). Our 
fitting method could deal successfully only with in- 
stances of sufficiently separated spectral features, 
where the local minima drop below 50% of the 
peak. More blended features are straddled with a 
single broad profile fit. From a careful examina- 
tion of individual cases of blended and unblended 
profiles, based on the presence or absence of in- 
tersecting linear position-velocity features as de- 
scribed above, we established a cut-off in fit qual- 
ity of = 25 that most effectively seperated these 
two regimes. Only those fits with a below this 
cut-off were retained, resulting in rejection of some 
4% of the total. Concentrations of rejected fits are 
located in the expected regions (noted above). 

The results of our spectral fitting were recorded 
as images of the physical parameters, N^f, Ts and 
a NT together with the integral of brightness tem- 



perature over velocity that corresponds to the col- 
umn density of the fit, / Tp*dV. This allowed cal- 
culation of a total corrected column density image 
from, 



Nj^f = iV^f +1.823x10^^ 



T'E°^dV- 



Tp'dV 



(5) 

where / T'^°^dV is just the usual image of inte- 
grated emission. In this way any spectral com- 
ponent not accounted for by the fit to the peak 
was retained, albeit with the assumption of neg- 
ligible opacity for this residual. Those lines-of- 
sight with insufficient fit quality (x^ > 25) were 
assumed to have negligible opacity. The two col- 
umn density images (assuming negligible optical 
depth and fitting for it) are contrasted in Figs. fTSl 
and [m Substantially more fine-scale structure is 
apparent in the opacity-corrected column density, 
primarily in the form of compact features and nar- 
row filaments. The narrow filaments, in particular, 
correspond to those seen in the images of velocity 
coherence, shown in Figs.[6]-[8l Peak column densi- 
ties are enhanced by about an order of magnitude 
over the uncorrected version. 

Given the uncertainties in the fit parameters, 
we do not attach undue significance to each pixel 
in the images of corrected column density, since 
they may easily be in error by a factor of 2 or 
more, but are encouraged by the spatial continuity 
of the solutions, since the fit to every line-of-sight 
is carried out independently of all others. Some 
indication for the relevance of the fit results can 
be obtained by comparison with high signal-to- 
noise observations of H I absorption toward back- 
ground continuum sources. As noted previously, 
several bright backgro und cont inuum sources that 
have been studied by iBraun fc Walte rboa (1992) 
fall fortuitously near regions of significant self- 
opacity, particularly the sources B0039+412 and 
B0044+419 for which peak optical depths of 
Tmax = 0.46 ± 0.02 and Tmax = 1-06 ± 0.06 are 
observed. (Note that (1 — e"'^) rather than r 
is tabulated in the reference.) The model fits 
show very substantial structure in the vicinity of 
these background sources, making it difficult to 
estimate the likely line-of-sight parameters. We 
have determined the average peak model opacity 
over a 3x3 pixel region that is offset azimuthally 
from the direction of each of these sources by one 
beamwidth. For these two lines-of-sight this yields 
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Tmax = 0.4 ± 0.4 and r^ax = 1-1 ± 0.8. Although 
the pixel-to-pixel variation is very large, and may 
well represent real structure rather than fitting 
errors, the order-of-magnitude of the opacity is 
estimated correctly in these cases. 

The total detected H I mass of M31 is increased 
by the opacity corrections by about 30%, from 
5.76x10^ to 7.33 xlO^ Mq. In view of the neces- 
sity of discarding all blended profiles, this should 
be regarded as a lower limit to the actual mass cor- 
rection (within the limitations of our approach) . 

The best-fit spin temperature and non-thermal 
velocity dispersion are shown in Figs. [T7HT51 Re- 
gions of significant opacity {\og{NHi) > 22) have 
spin temperatures of less than about 100 K and 
are often organized into patchy filamentary com- 
plexes. While there is some organized spatial 
structure in the distribution of non-thermal veloc- 
ity dispersion, it is not obviously correlated with 
either the total column density or spin tempera- 
ture. Not surprisingly, there is a significant corre- 
lation of the non-thermal velocity dispersion of the 
fits with the image of velocity coherence (shown in 
the right hand panels of Figs. [6H8]). 

4.3. Radial Trends 

The basic radial trends of observed and derived 
quantities are illustrated in Fig. [191 The peak ob- 
served brightness temperature (Panel A) is seen 
to first increase dramatically with radius out to 
about 12 kpc and then subsequently to decline 
out to the the largest detected radii. Even the 
median peak brightness exceeds 50 K near 12 kpc, 
while the distribution extends out beyond 100 K. 
This can be co mpared with th e very general trend 
documented by [Braun (1997) for a positive gra- 
dient in peak brightness to occur within the ac- 
tively star forming disk of nearby galaxies. The 
unprecented physical sensitivity obtained for M31 
allows us to track the atomic gas properties out 
to more than three times larger radii. The ap- 
parent column density distribution in Panel B fol- 
lows a similar trend, illustrating the exponential 
character of both the increase to- and the sub- 
sequent decline from- the peak near 12 kpc. A 
6 kpc scale-length exponential function is over- 
laid on the plots in Panels B, C, D and H to il- 
lustrate how well this scale-length describes the 
gas mass density. Intriguingly, the old stellar 
population of M31, as traced by 3.6 /im Spitzer 



imaging between 1 and 12 kpc radius is best- 
fit with a 6.08ib0.09 kpc scale-length exponential 
( Barmbv et al. 2006f l. The corrected H I column 
density, shown in Panel C, has a median ridge- 
line that closely tracks the apparent column, but 
also demonstrates that outliers to higher corrected 
column occur primarily at radii between about 7 
and 30 kpc. Adding in the molecular contribution 
to mass surface density (as described in detail in 
§ 14. 4p in Panel D yields only minor changes to the 
radial trend of the atomic column density. The 
most significant is a slight addition of mass den- 
sity near 12 kpc, where the molecular contribution 
permits the total to better track the 6 kpc expo- 
nential to slightly smaller radii. Inside of about 
12 kpc is a sharp cut-off that can be described 
with a 3 kpc exponential scale-length. Beyond 
about 28 kpc, where opaque H lis no longer found, 
there is a significant truncation of the disk expo- 
nential (neutral) mass density, corresponding to 
a 2 kpc scale-length. The corresponding face-on 
column density (scaled by l/cos(78°)) where this 
outer truncation sets in, is 5xl0^^cm~^. 

The radial distribution of velocity coherence 
(Panel E) shows a similar trend, rising rapidly 
out to about 10 kpc, then declining out to about 
30 kpc and then plummeting. Both the derived 
non-thermal velocity dispersion (Panel F) and spin 
temperature (Panel G) show some broad system- 
atic decline with radius beyond 10 kpc. For the 
spin temperature, it is only the lower envelope 
of the distribution which represents "opaque" fits. 
These "opaque" cool component temperatures are 
as low as 20 K near 7 and 25 kpc, but increase to 
about 60 K near 10 kpc. This is quite distinct from 
the most highly populated ridge-line, which repre- 
sents predominantly "transparent" lines-of-sight. 
The ridge-line shows a similar trend, but is offset 
to higher temperatures by about an order of mag- 
nitude. Current massive SFR density (as defined 
below in § l4.4p is shown in Panel H. The apparent 
nuclear peak in SFR density is due to other pro- 
cesses (as discussed below), while current coverage 
only extends to about 20 kpc. Current SFR den- 
sity is highly peaked near 12 kpc, where gas mass 
densities are highest and does not track the 6 kpc 
exponential of the median gass mass density. 
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Fig. 15. — Column density of H I assuming negligible opacity of the central 25% of the survey region. 
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Fig. 16. — Column density of H I derived from a fit to each spectrum of the isothermal model described in 
]for the central 25% of the survey region. 
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Fig. 17. — The logarithm of the spin temperature of H I derived from a fit to each spectrum of the isothermal 
model described in HA.2\ for the central 25% of the survey region. 
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Fig. 18. — The logarithm of the non-thermal velocity dispersion of H I derived from a fit to each spectrum 
of the isothermal model described in ^^4.21 for the central 25% of the survey region. 
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4.4. The Surface Density of Gas Mass and 
Star Formation Rate 

A particular form of correlation of observable 
parameters deserves special consideration, namely 
that between the surface density of gas mass 
and the star formation rate. Many authors have 
addressed the relationship between these prop- 
erties on both a global and loc al level, begin- 
ning with the pioneer ing study of ISchmidd (l959l) 
and most rcccnth r bv lKennicutt et all (|2007h and 
Irhilker et al] (2007h . who probe the resolved re- 
lationship of these parameters down to scales of 
500 pc within NGC 5194 (M51a) and 400 pc within 
NGC 7331 respectively. Given the proximity of 
M31 we can extend this work by more than an 
order of magnitude to smaller surface area at spa- 
tial scales of 100 pc, albeit with a substantial 
inclination correction that introduces some poten- 
tial linear smearing (by a factor of 4.8) in one 
dimension. Mass and star formation rates have 
been calculated after background subtraction and 
spatial smoothing of all relevant images to the 
same 30 arcsec (113 pc) Gaussian resolution as 
applies to the H I properties. 

The mass surface density was calculated from 
the integrated CO(l-O) emission obs e rved with 
the IRAM telescope by iNieten et all (|2006t ) as- 
suming a mean atomic mass of 1.36 uyh and an 
Nffa/Ico conversion factor of 2.8x lG^°cm^^/K- 
km s~^. This differ s from the procedure adopted 
by iKennicutt et all (|2007t ). (hereafter K07) who 
keep the comparison in terms of hydrogen nuclei 
(rather than including a Helium mass correction) . 
The helium mass correction amounts to a shift by 
-1-0.13 dex in mass density. Note that the CO 
data only extend over the central 2x0.5 degrees 
(or 28x32 deprojected kpc in diameter). Radii 
beyond about 15-20 kpc therefore only have an 
H I contribution to the calculated mass density. 
Although the CO brightness at the edges of the 
sampled region is generally faint and has a low sur- 
face covering factor, isolated pockets of CO may 
well be present at larger radii. 

The star formation densi ty was calcula, t ed fol- 



lowing the methodology of iThilker et al.l pOOTl ) 
from the Spitzer I RAC 8/im and MIPS 24jnm sur- 



face brightness of iGordon et al.l (120061) together 



metric" luminosity was assumed to be, L[bol) = 
vpuv Lu.,obs{FUV) + {l — r])L{IR) after correcting 
the FUV image for foreground extinction by the 
Galaxy, and assuming 77=0.32. The term "bolo- 
metric" in this context is a misnomer, since the 
quantity is not assumed to represent a true inte- 
gral of the entire SED, but merely to provide sensi- 
tivity to both directly emitted UV and re-radiated 
IR light. The value of ry, which represents the frac- 
tion of the IR luminosity due to an older stellar 
population, is expected to vary substantially with 
environment within a galaxy. Ideally, this would 
also be modeled in a position dependent way, but 
that is beyond the scope of the current discussion. 

The total (3-llOOum) IR luminosity per pixel 
was estimated using the IRAC/MIPS data and the 
relation, 



log 



L{IR) 



1.06 + 0.475 log 



with the GALEX FUV im a ge of iThilker et al 
(l200,5[ ) and iGil de Paz et all (|2n07h . The "bolo- 



(6) 

and then converted to a star formation rate den- 
sity using the projected pixel area and the SFR 
calibration: 

log SFR{bol){MQyr-^) = log L{bol){LQ) - 9.75. 

(7) 

The 24(1111 data only extend over the brighter disk 
region of ^3x1 degrees. 

Only pixels with significant emission in the spa- 
tially smoothed images were used for the compar- 
ison, which resulted in thresholds at a brightness 
of 0.7 K-km s~^ for the CO emission, H I columns 
greater than 19 in the log and sky-subtracted IR 
brightness greater than 0.06 MJy/sr at 24/im and 
0.1 MJy/sr at 8/im. As can be seen in Fig. [T51 
the spatial coverage and sensitivity in the IR only 
extends out to a radius of about 20 kpc. 

We contrast three different versions of the star 
formation law in Fig. [20] utilizing the molecular 
only, the molecular plus apparent H I and the 
molecular plus opacity-corrected H I column den- 
sity. The molecular-only plots are overplotted 
both the relation of K07 wit h slope 1.37, a value 
which has also been found bv lHever et al. ( 2004^ 
in M33, as well as a line with unit slope corre- 
sponding to a c onstant gas de p letion time-scale, 
as advocated by Wong fc Blit d ( 20021 ). (hereafter 
WB02). The slope of the molecular-only relation- 
ship is approximately matched by the simple lin- 
ear dependency. The gas depletion time-scale of 
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the plotted line is 1 Gyr (for a N/f^/Ico conver- 
sion factor of 2x 10^°cm~^/K-km s~^ as assumed 
by Wong & Blitz). This is similar to many of 
the galaxies analyzed in WB02, although occur- 
ring at mass and star formation densities that are 
at least an order of magnitude lower within M31. 
The M51a molecular-only relation of K07 is both 
steeper and shifted to higher mass densities by at 
least 0.3 dex. 

The total-gas mass density plots are over- 
plotted with the relation fit by K07 to the 
M51a data with slope 1.56. Note that we have 
plotted star formation rate density in units of 
Mopc^^Gyr^-^ rather than the Mokpc^^yr^^ of 
K07 implying a shift of 3 dex. The different units 
of both axes (relative to K07) have been taken 
into account for the overplotted curves. Remark- 
ably, the M51a relation provides a very good fit 
to the M31 data, despite the fact that it has been 
derived in the strongly molecular-dominated inner 
disk of M51 and is being applied thoughout the 
atomic-dominated disk of M31. The highest con- 
tours of correlation probability are about twice as 
narrow for the total-gas relation compared to the 
molecular-only one. In the relation using apparent 
H I, there is a very steep truncation to mass den- 
sities exceeding about 10 M0pc~^, together with 
a saturation that sets in near 16 Mopc"^. This 
truncation disappears in the opacity-corrected H I 
relation, together with a small shift of the peak 
correlation ridge-line to higher mass densities. 
The saturation effect near 16 Mopc^^, is also 
alleviated although not entirely eliminated. 

Data at the smallest radii (<5 kpc) yield the 
diffuse region of elevated apparent star formation 
in the total gas panels of the figure. The enhanced 
24/im emission from the nuclear regions of M31 
is accompanied by both recombination- and for- 
bidden emiss i on lin e gas (as first documented by 
Jacobv et all (|l985h ) with line ratios suggestive of 
shock heating, and as such is not due to massive 
star formation. Essentially no CO is detected from 
the nuclear region, so this feature is not apparent 
in Panels A and D. Beyond about 5 kpc, the same 
basic correlation is seen to apply at all radii, with 
only the range of the observables showing some 
systematic variation, but not the basic ridge-line 
of their correlation. The sharp truncation of the 
molecular-only mass density near 1 Mqpc^^ is a 
consequence of the limited sensitivity. The same is 



true for star formation rate densities below about 
0.2 Mgpc-^Gyr-i. 

At star formation rate densities below about 
0.4 Mopc~^Gyr~^ and total gas surface densi- 
ties of about 5 Mqpc~^ there appears to be a 
departure from the power-law relationship. The 
same feature can be seen in the NGC 7331 data 



of Thilk er et al.l (|2007r ) at a similar location, al- 
though this is near the noise floor in both cases. 
This may indeed be a distinct threshold for the for- 
mation of molecular clouds from which stars will 
form. The most substantial contribution to this 
possible threshold comes from radii in excess of 
15 kpc. 

Perhaps the most surprising aspect of Fig. is 
the good statistical correlation of total gas mass 
and star formation density down to spatial scales 
of about 100 pc and gas masses of only STO** M©, 
well inside the regime of individual, giant molecu- 
lar clouds. 

5. Conclusions 

Our high resolution and signal-to-noise study of 
the neutral ISM of M31 has permitted a number 
of novel results: 

1. We detect ubiquitous self-opaque H I fea- 
tures, discernible in the first instance as fila- 
mentary local minima in images of the peak 
H I brightness temperature. Such local min- 
ima are often accompanied by systemati- 
cally broader line profiles, apparent in im- 
ages of the velocity coherence. High opac- 
ity features are organized into complexes of 
more than kpc length and are particularly 
associated with the leading edge of spiral 
arms. We suggest that such features are the 
counterpart of large-scale H I self-absorption 
(HISA) features seen in the Galaxy. They 
are the consequence of a "sandwich" geom- 
etry of two warm, semi-opaque layers which 
flank a colder, opaque core. Such a geome- 
try leads to widespread HISA in the edge-on 
Galactic disk, but is spatially resolved into 
distinct parallel filaments in the inclined disk 
of M31. Just as in the Galaxy, there is only 
patchy correspondence of self-opaque/HISA 
features with CO(l-O) emission. 

2. We have produced images of the best- 
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Fig. 19. The radial distribution of several diagnostie properties of M31. The plotted quantities are the 
peak H I brightness temperature, apparent H I column density, corrected H I column density, the face-on 
gas mass surface density in M0pc~^, H I velocity coherence, the best-fit non-thermal velocity dispersion 
(labelled log(s_NT)) the best-fit cool H I temperature (labelled log(T_S)) and the estimated face-on star 
formation rate density in M0pc~^Gyr~^ based on FIR and UV surface brightness. Note that the central 
peak in panel H is not due to star formation (see text). Contours are drawn at factors of two in correlation 
density beginning at 0.1% and ending at 51%. The greyscale extends from and 100% with a square-root 
transfer function. The gray solid line in panels B, C, D and H is an exponential with 6.0 kpc scale-length. 
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Fig. 20. — Star formation rate density versus gas mass density in M31. The face-on gas mass surface density 
has units of Mqpc~^ and the face-on star formation rate density, MQpc^'^Gyr^^. The mass surface density 
is calculated with molecular only-, molecular plus apparent H I- and molecular plus opacity-corrected H I 
column density from left to right. The data at all radii are shown on the top, while the radially binned data 
is shown on the bottom. Contours are drawn at factors of two in correlation density beginning at 0.5% and 
ending at 64%. The greyscale extends from and 100% with a square-root transfer function. The colors red, 
green and blue correspond to bins in galactocentric radius of 0-8, 8-16 and 16-24 kpc. The diagonal lines of 



slope 1.37 (molecular only) and slope 1.56 (molecular plus atomic) are the relations given in lKennicutt et al 
(|2007l ). A second line of unit slope is plotted in the molecular-only plots, corresponding to a gas depletion 
timescale of 1 Gyr. 
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fit physical parameters; spin temperature, 
opacity-corrected column density and non- 
thermal velocity dispersion, for the brightest 
spectral feature along each line-of-sight in 
the M31 disk. This represents a major step 
forward in the quantitative assessment of H I 
opacity in a galactic context. Opaque atomic 
gas is organized into filamentary complexes 
and isolated clouds down to the resolution 
limit of 100 pc. The spin temperature of 
opaque regions first increases systematically 
with radius from about 20 to 60 K (for radii 
of 5 to 12 kpc) and then smoothly declines 
again to 20 K by about 25 kpc. Opacity 
corrections to the column density exceed an 
order of magnitude in many cases and add 
globally to a 30% increase in the atomic gas 
mass over that inferred from the integrated 
brightness under the usual assumption of 
negligible self-opacity. It will be particularly 
interesting to undertake a similar analysis of 
more nearly face-on systems like M33 and 
the LMC, where spectral blending will be 
even less of an issue. 

3. We have constructed the radial distribution 
of gas mass density which illustrates a re- 
markably extended exponential decline with 
6 kpc scale-length between 12 and 28 kpc. 
Inside of 12 kpc there is a rapid decline in 
neutral mass density with 3 kpc scale-length, 
while beyond 28 kpc at a face-on column 
density of 5 x lO-'^^cm"^ is a truncation of the 
atomic disk with 2 kpc scale-length. Perhaps 
fortuitously, the extended disk scale-length 
of 6 kpc matches that of the old stellar pop- 
ulation traced by 3.6 fim emission. 

4. We have extended the resolved correlation of 
star-formation-rate- with gas-mass-density 
down to the smallest physical scales yet 
reached, of only 100 pc and lO'* Mq; more 
than an order of magnitude smaller in area 
and mass than has been possible previously. 
Unlike other galaxies, in which the gas mass 
is dominated by molecular gas as traced by 
CO emission at small radii, M31 is domi- 
nated by atomic gas at all radii. The re- 
lation of molecular-mass- to star-formation- 
density has a large dispersion, but has a 
slope near unity and a corresponding gas 



depletion timescale of about 1 Gyr, similar 
to ma ny galaxies studied by IWong fc Blitz 
((2002) • The relation between total-gas- 
mass- and star-formation-rate-density is sig- 
nificantly tighter and is fully consistent 
in both slope and normalization with the 
relation found i n the molecule- d omina ted 
disk of M51 by iKennicutt et al.l (|2007f ) at 
500 pc resolution. Together, the M31 and 
M51 data demonstrate the same relation- 
ship spanning mass-densities from about 5 
to 700 Mqpc"^. Use of opacity-corrected H I 
columns yields a more symmetric distribu- 
tion with less evidence of truncation effects 
than with the apparent H I column. Below 
about 5 Mqpc~^, there is a down-turn in 
star-formation-density which may represent 
a real local threshold for massive star for- 
mation. The corresponding threshold cloud 
mass IS about S-IO* Mq. 
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